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Abstract 

We consider in the soft-collinear effective theory semi-inclusive hadronic B decays, B — > XM, 

in which an energetic light meson M near the endpoint recoils against an inclusive jet X. We 

focus on a subset of decays where the spectator quark from the B meson ends up in the jet. The 

branching ratios and direct CP asymmetries are computed to next-to-leading order accuracy in 

a s and to leading order in l/m?,. The contribution of charming penguins is extensively discussed, 

and a method to extract it in semi-inclusive decays is suggested. Subleading l/m& corrections and 

^ ' SU(3) breaking effects are discussed. 
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I. INTRODUCTION 



ess attention over the 



Semi-inclusive hadronic decays B — > XM have received much 
years in contrast to the widely studied exclusive two-body B decays 
As we will show in this paper, semi-inclusive hadronic B decays in the endpoint region, 
where M is an isolated energetic meson and X is a collinear jet of hadrons in the opposite 
direction, are theoretically simpler than the exclusive two-body B decays in many respects, 
yet still address many of the questions that had been debated in the context of the two-body 
B decays. Using the soft-collinear effective theory (SCET) [ill EH I2-E3 predictions of 
semi-inclusive decays can be improved systematically and lead to the following advantages. 
Firstly, larger data samples can be included by considering inclusive jets with a variety of 
final-state particles forming the collinear jet. Secondly, as in exclusive B decays [14, 15]. 
the four-quark operators in the weak Hamiltonian factorize at leading order in l/m^ into 
a product of a heavy-to-light current and a collinear current, with no strong interactions 
between these two currents. Thirdly, the inclusive collinear jet is described by the jet 
function, which is obtained by matching the full theory onto SCETj at the scale p 2 x ~ m^A 
with A ~ Aq CD . Since the same jet function appears at leading order in B — > X s ^ or 
B — > XlV inclusive decays, many of the hadronic uncertainties cancel by taking ratios. 
Finally, the contribution of charming penguins can be implemented systematically using the 
effective theory. Studying B — > XM decays can thus offer a theoretical handle to probe 
nonperturbative effects of charming penguins. 

In order to see these advantages clearly, we consider the decays B — > XM in which 
the spectator quark of the B meson goes to the inclusive jet. It is straightforward to treat 
other decay modes without this constraint 16], but would involve more calculation including 
spectator interactions, and we do not discuss it further here. The decay rate for B — > XM 
at leading order in l/m?, can then be schematically written as 

<,[ {\T®<p M \ 2 )-{J®f) } (1) 



dE 



M 



where T is a collection of hard coefficients obtained in matching full QCD onto SCETj and 
J is the discontinuity of the jet function describing the fluctuations of order m b A in the 
forward scattering amplitude of the heavy-to-light currents. <pM and / are the light-cone 
distribution amplitude (LCDA) for the light meson M and the .B-meson shape function, 
respectively. The ® sign implies the appropriate convolution. The convolution J ® / 
in Eq. is universal in the sense that the same convolution appears in B — > X s 7 and 
B — > X u lV decays n], 17]. Therefore, if we take the ratios of the decay rates for B — > XM 



and, say, the decay rate for B — ► A s 7, this convolution cancels out and the only surviving 
nonperturbative parameters are the LCDAs. 

Another interesting but complicated problem common to two-body B decays and B — > 
XM decays in the endpoint region is the contribution of intermediate charming penguins, 



which can be of nonperturbative nature [18|. There has been a disagreement on how to treat 
this contribution between the recent SCET analysis of the two-body B decays 0, E3] and 
the approach of QCD factorization j^jj. The question is whether or not the long-distance 
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effects of charming penguins are of leading order in l/m&. Long-distance contributions arise 
when intermediate charm quarks lie in the non-relativistic QCD (NRQCD) regime with 
small relative velocity v*. These contributions are of the form a s (2m c ) f (2m c / m^v * 2~3 |. 
where /(2m c /mj) is a factor which accounts for the phase space in which the charm quarks 
have small relative velocities. In QCD factorization 0, the claim is that the phase 
space suppression near the threshold region is strong enough so that the nonperturbative 
contributions are subleading. On the other hand, Bauer et al. 1*2 *2f ] argue that since 
2m c /mf, is of order unity so is f(2m c /mb), and the nonperturbative contribution of charming 
penguins can be of leading order. In this paper, we suggest how to resolve the issue of 
charming penguins in B — > XM decays. If the nonperturbative contributions of charming 
penguins are really suppressed, then the decay rates at leading order in 1/mb are completely 
determined in terms of the perturbatively calculable hard kernels convoluted with LCDA, 
once normalized to the B — > X s 7 rate. If nonperturbative charming penguins are not 
suppressed, they will show up experimentally as a sizable deviation from purely perturbative 
predictions, which we will discuss in detail. 

The paper is organized as follows: In Section |H] we describe the kinematics for B — > MX 
decays. In Section IIII[ we set up the operator basis for the decays B — > MX and compute 
the radiative corrections at next-to-leading order (NLO) to derive the renormalization group 
equations for the operators. In Section IIVI we present a factorized form for the semi- 
inclusive B decays in the endpoint region. Section El is devoted to the contribution of 
charming penguins, considering two possible scenarios in which the charm quark is regarded 
as either hard-collinear or heavy. The contribution of charming penguins in the heavy quark 
limit mj, m c — > oo with m c /mb fixed is considered in detail. In Section EH we discuss the 
corrections to the leading order prediction, including SU(3) breaking effects. In Section 
IVII[ we perform the phenomenological analysis of B — > MX decays and predict the decay 
rates and CP asymmetries. The method to extract the effect of charming penguins is also 
discussed. We conclude in Section IVIIII In Appendix A we present the Wilson coefficients 
for the operators at NLO in SCETi. In Appendix B the detailed analysis of charming 
penguins in the heavy quark limit is discussed. 



II. KINEMATICS 

Using SCET, solid predictions can be made for hadronic semi- inclusive B — > MX decays 
in the endpoint region. In the rest frame of the B meson, the outgoing energetic meson 
M with p 2 M ~ A 2 moves in the rT M direction, while the inclusive hard-collinear jet with 
p 2 x ~ Ami, is in the direction, where n 2 = n 2 = 0, n ■ n = 2. We can choose the reference 
frame in which the transverse momentum of M is zero. The momenta p^ M and p\ can be 
written in terms of the light-cone coordinates p^ = (n ■ p,n ■ p,p±) as 

Pm = (o, n • p A f , o) +0(A 2 /m B ), 

Px = (m B ,m B - n ■ p M ,0) + 0(A 2 /m B ), (2) 
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with p B = tubv^ = p M + p x , where 2u M = n M + n M . We consider the endpoint region in 
which m B — n ■ pu ~ A, so that p\ ~ m B A. 

At the quark level, the b quark has momentum p 1 ^ = m b v fl + P, where / M is the residual 
momentum of order Aqcd- The b quark decays to a quark-antiquark pair moving in the 
n direction which hadronizes into the meson M, and another quark with momentum pj 
moving in the n direction, which combines with a spectator antiquark to form the outgoing 
jet X. The momentum pj can be written as (dropping terms of order A 2 /m b ) 

= m bV » + F- p M = m b — + m b (l - x M ) y + P = m b — + F, (3) 

where xu = n ■ PAi/rrib = lEulm^. In the endpoint region, 1 — xm ~ K/m b . Since the 
invariant mass squared p 2 j of the jet is time-like, the range of the residual momentum fc M in 
Pj is < n ■ k < n ■ px ■ Since the residual momentum of the heavy quark n ■ I is smaller 
than A = rriB — rn b , the region of n ■ I, which has support for the S-meson shape function, 
is 

-m 6 (l - x M ) < n ■ I < A. (4) 



III. MATCHING AND EVOLUTION IN SCETi 



We employ a two-step matching in computing and evolving the hard coefficients. First 
we construct the operators for the decays in SCETj by integrating out degrees of freedom 
of order m b . The Wilson coefficients of the operators are obtained by matching full QCD 
onto SCETi. The decay rates of the semi-inclusive B decays are obtained from the forward 
scattering amplitude of the time-ordered product of the heavy-to-light currents, as shown 
in Fig. [TJ In the next step, we match SCETi onto SCETn by integrating out the degrees of 
freedom with p 2 ~ m b A. As a result, the jet function is obtained, the discontinuity of which 
contributes to the semi-inclusive hadronic B decay rates. 

The effective weak Hamiltonian in full QCD for hadronic B decays is given as 



w 
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FIG. 1: Tree- level diagram for the forward scattering of the heavy-to-light currents in SCETi whose 
discontinuity gives the semi-inclusive hadronic B decay rates in the endpoint region. The double 
lines denote a heavy quark, the intermediate line is the hard-collinear quark in the n direction with 
m 2 (l — xm), while the upward collinear quarks move in the n direction. 
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where the operators are 

0\ = (pb) V - A (qp)v- A , O p 2 = (j)pb a )v- A {q a Pp)v-A, (6) 

C 3 ,5 = (qb)v-A ^(q'q')vTA, C 4 ,6 = (qpb a ) v _ A ^(7 a g^)y T A, 

q> q> 

7 ,9 = (qb)v-A -^-(q'q')v±A, 08,10 = (qpb a ) V - A ^ -ifiValp^A, 

q> q> 

Oy = + 7 5 )&, o g = -S^g;t°(i + 75 )6. 



Here = V p bV* q is the CKM factor and V ± A = 7 M (1 ± 75). The summation over g' 
includes u, d, s, c and b quarks. Operators with q = d (q = s) describe the AS = (AS = 1) 
effective weak Hamiltonian. 

The effective Hamiltonian in SCETi at leading order (LO) in f/m& is (with charm quarks 
integrated out, nonperturbative charm contributions will be discussed in Section IV1) 

^^e^ec?®^, (7) 



quarks 



p=u,c i=l 



where the relevant four-quark operators are 

C 1 = [u n fP L Y^b v ] [q n $P L Un] u , C 2 ,3 = [q n ^P L Y^b v ] [Un4 P L,RUn] u , (8) 

The summation over </ includes w, <i and s quarks and Pl,r = (1 =F 75)/2. In Eq. (J7J), ® 
denotes the convolution 

Cf®0i= / duC?(u)0i(u), (9) 







and the subscript u in Eq. (jHJ) refers to the variable in a 5 function, which is defined as 

r / n • 1 

[qSPLOn] u = [Qn \U - "2^- ) i P L(ln\ ■ (10) 

The q n and are the gauge-invariant quark fields 

Qn = Wi&>, qn = Wl&\ (11) 

given in terms of the collinear fermion fields £4 , ^i 9 "* of flavor q and the collinear Wilson 
lines W n , W n in the n and rT directions, respectively. The ultrasoft (usoft) Wilson line in 
the n direction, Y n , is obtained after redefining the collinear fields to decouple collinear and 
usoft degrees of freedom [l3j | . 

There are also color-octet operators corresponding to the operators in Eq. (J7|), e.g., 

O x (u) = [u n Y^Y w fP L T a Yib v ] [q^P L T a UjT ] u , (12) 



but the matrix elements of the octet operators between hadronic states vanish and are 
therefore not relevant here. The Wilson coefficients Cf(u) in Eq. (J7J) encode physics at the 
hard scale my, and are perturbatively calculable in powers of a s {m^). They are known at 
NLO in a. 



la, |2l( and are listed in Appendix |XJ Note that the Wilson coefficients Cf(u) 



exhibit nonzero strong phases at NLO from integrating out the intermediate on-shell quarks. 

In matching SCETi onto SCETn at /i ~ y/rribA, the operators in the Hamiltonian in 
Eq. (J7|) are first evolved down from rrif, to yU using the renormalization group (RG) equation 
in SCETi. The operators in Eq. (JZj) factor into a heavy-to-light current and a collinear 



current as 



0(u,n)= q n T H Y^b v q„T C 9 



(13) 



where Yh,c are the Dirac structure in each current. There are no strong interactions between 
the two currents to all orders in a s in SCETi. At order a s , the radiative corrections in Fig- EI 
show explicitly that this is true. As a result, the operator O is multiplicatively renormalized, 
and there is no mixing between color-singlet and color-octet operators due to factorization. 
The renormalized operator Or and the bare operator Ob are related by 

O r (u,h) = J dvZ' l (u,v ,h)Ob(v,h) = J dvZ^^Z^fav, /jl)O b (v, fi), (14) 

where the counterterm Z is a product of the counterterms Z^ and Zc from the radiative 



b) 



c) 




FIG. 2: Feynman diagrams at order a s for the four-quark operators in SCETi. Note that there is 
no strong interaction between the two currents. 



04 is a sum over a product of currents. When considering the spectator quark going into the jet, only 
one of the terms will contribute. 
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corrections of and Jq. This leads to the RG equation 

= -1hJh(h)Jc{u,ijl) - Jh{h) J dvy c (u,v, n)J c (v, fx) 
whereas the RG equation for the Wilson coefficients is written as 



(15) 



djji 



dv 



(16) 



At next-to-leading logarithm (NLL), the anomalous dimensions for Jh and Jc are given by 



>f . a a s Cp /5 a 
— In 



71 



\2 n 



m b 



a s C F 



2tt 



5(u — v) + 2— ( 1 + — ^ 9(v — u) + ( u, v <-> u, v) 

v V (v — u) + / V / 



(17) 



, (18) 



where Cp = (N 2 — 1)/(2N) with N the number of colors. The subscript '+' denotes the 
plus distribution, and u = 1 — u. The one-loop result for in SCET was first obtained 
in Ref. [loj . while the part of the two-loop result containing a 2 s ln(///m&) needed at NLL 
accuracy has not yet been calculated in SCET. Extracting it from the full QCD calculation 

avors. The one-loop 



24j . one gets B = iV(67/18 — vr 2 /6) — 5n//9, where n/ is the number of f 

ull QCD hard kernel calculations 



result for 7c can be taken from the 
with the determination in SCET [27 1 
At one loop, Eq. (fTo^) can be written as 



251 126| , which agree 



djj, 



2vr 



a s Cp 1 
2n uu J 



dvV(u,v)Ci(v, fi), (19) 



where V(u, v) is the Brodsky-Lepage kernel |2 



V(u,v) 



uv (l + - — ^ 9(u — v) + (u, v u, v) 

v [u — v) + y v / 



(20) 



The eigenfunctions of Eq. ()19j) are given by the Gegenbauer polynomials Cn 2 {2u — 1), which 
satisfy 

/•i 

cW(w, f )Cf 2 (2w - 1) = e n wwCf 2 (2w - 1), (21) 



with the eigenvalues 



-4 + 



(n + l)(n + 2) 



n+l ^ 

4 ^P 



(22) 



fc=2 



since the four-quark operators in SCETj are partially governed by the light-cone conformal 
symmetry with the highest weight of the conformal spin j = 2 + n ^2J^ . 
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We can now expand the Wilson coefficients in terms of the Gegenbauer polynomials, 

C l (u,v) = Y,C 3 n / \2u-l)Anri, (23) 



n=0 



a virtue of which is that A™ with different n do not mix to one loop. The solution of the 
RG equations for A™, 



d A n( \ \ i \ Q a sCF oc s C F 



yield the Wilson coefficient at the scale /i to NLL order 

b{u,n) = ^Cf 2 (2 M - 1) A?(m b )exp[In(j*,m b )], 



(24) 



(25) 



with 



I n (n,m b ) 



An C F 



1 In z 

z 

2B 



a s {m b ) (31 

C F (\\ \, -Hi, 



+ pI Cf 



— 1 + z — In z + — In 2 z 



1 — z + In z 



(26) 



where z = a s (fi)/a s (m b ). The coefficients of the QCD (3 function are (3q = (11 — 2n/)/3, and 
Pi = 34N 2 /3 — 10Nrif/3 — 2C F nj. From the orthogonality of the Gegenbauer polynomials 
the coefficients A™(m b ) are 



A?(m b ) 



4(2n + 3) 
(n + l)(n + 2) 



/ du u(l -u)Ci{u,m b )Cl /2 {2u- 1). (27) 
Jo 



At NLL, only LO values of the Wilson coefficients at // = m b are needed. Since these are 
independent of the momentum fraction u, we have A™{m b ) = Ci,Lo(^b)<W 



IV. SEMI-INCLUSIVE DECAY RATES 

The decay amplitudes for B — > XM, in which a spectator quark of the B meson ends up 
in the jet X, are schematically 

(XM\Hj\B) = ^ J\uT^\u^) (M\[^P L q£jO)(X\qjkP L YX\B), (28) 

where the hard kernel T$ is given by the sum of the products of the CKM factors \$ = 
V pb V* and the Wilson coefficients Cf. Here q denotes the flavor of the outgoing quark in the 
heavy-to-light current. The matrix elements for the meson M are related to the LCDA by 

(P|[^75^]J0) = -2tf P E P Mu^), (29) 
041 WM]J0) = 2if v EvMu,»), (30) 
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Mode (LAD = 1J 




Mode [flj = UJ 


1 M 


If + 


Ap (L l + L 4 J 


<n- — V + V + 
71 A / jj YY 


Ap \yi -f- L- 4 ) 


7T { * )0 x- 






"p 1-4 




A p L A — A t (L.5 (-6 J 




— A t + 1-6 J 




^(A p SJ Cf + X { t s) C 3 ) 


n°X° s 




P°X% 


^(A?Cf-A^C 3 ) 


(PX° 


-^(Ap L 2 - A t L 3 - Ap L 4 J 


"X% 


^-Cf-^-(C 3 + 2C 5 + 2C 6 


) uXl 


.(d) .(d) 

^-(Cf + Cf)-^-(C 3 + 2C 5 + 2C 6 ) 




A«(C? + Cf) 


ir-Xf/p-Xf 


Ap^ + Cf) 



TABLE I: Hard kernels T^' for 5 /5~ -► XM (above horizontal line) and B s -> IM decays 
(below horizontal line), where only the strangeness content of the inclusive jet is shown. The 
summation over p = u, c is implied. The NLO Wilson coefficients Cf are given in Appendix lAl 



where P and Vl denote pseudo-scalar and longitudinal vector mesons respectively. Trans- 
versely polarized mesons, Vp, do not contribute at leading order, as in exclusive two-body 
B decays (for charming penguins, see below). Thus the decay amplitude can be written as 

(XM\H I \B)=i^f M E M [ duT$\u,fJ,)<p M (u,li)(X\qJlP L YX\B), (31) 

Jo 

where the hard kernels T^f for various processes are listed in Table 
In order to obtain the decay rates for B — > XM 



dT 



dE 



M 



(27T) 



m B 



dur! q) 



M VM 



Y,\(X\lnfPLYX\B)\ 2 &\PB-pX 



Vm) 



(32) 



x 



we use the optical theorem to relate the decay rate to the imaginary part of the forward 
scattering amplitude. We therefore consider the time-ordered product of the heavy-to-light 
currents 

T(E M ) = — [d 4 z e- ipM - z (B\T4(z)J H (0)\B), (33) 

TTT'B J 

where Jh( z ) — e l ^~ mbV ^' z q r $P L Y 1 lb v (z). Since there are no collinear particles in the B 
meson, the time-ordered product of the collinear fields can be written as 



(Q\Tq n (z)-q n (Q)\0)=ip(n-z)5 2 (z ± : 



dn ■ k 



- in - kn ^ 2 J P (n- k + ie), 



(34) 



which defines the jet function Jp = Jp(n ■ k) with the label momentum P. In SCETn, the 
remaining matrix elements are parameterized in terms of the B meson shape function, 



/(n-0 = 2 



dn ■ z 

4-7T 



„—in-ln-z 12 / tj \ 

e {±} v \ 



b Y 



(n-z/2) 



Y ] b 



(o)|£„) 



■(B v \b v Y n 8(n-l-n-id)YX\B v ), 



(35) 
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and the time-ordered product in Eq. can be written as 



T(E M ) = -2 j dn-l f(n- 1) J P (m b (l - x M ) + n ■ I + ie\ (36) 

J-m b (l-x M ) 

with the limits on n ■ I included according to Eq. (J3J). Taking the discontinuity, we obtain 
—IyrT^Em) — 2 / dn ■ I f \n ■ I) lmJp(m b (l — xm) + n ■ I + ie) (37) 

J-m b (l-x u ) /J 

m b 

where the nonperturbative function S is defined as the convolution of the B meson shape 
function and the imaginary part of the jet function. 

Combining Eqs. (jSH) and (jSZj), the factorized differential decay rate for the B — > XM is 

iB^XM) = -^f 2 M m 2 b x 3 M S(x M ,fio)H^(m b ,fi ), (38) 



dE M 8vr 



where is the convolution of the hard kernel and the LCDA, 



pi 2 

H^(m b ,^)= / duT$\u, /j, )(I)m(u, Ho) ■ (39) 
Jo 



The information on the LCDA can in principle be extracted from experimental data on other 
hard processes, while H$ can be computed in perturbation theory. It is worth mentioning 
that Eq. (j3~Hj) is independent of /z and fi. S(xm,Ho) is the convolution of the imaginary 
part of the jet function, which is computed in matching between SCETi and SCETn at fi 
and evolves down to fi, with the 5-meson shape function, evaluated at fi. The dependence 
on the low scale \i cancels between the two. In , T^f evolves from m b to /x 0) an d the 
LCDA (pM, which is the matrix element of the collinear quark operators, are evaluated at 
/i . The dependence on /i in will then cancel against ImJp. Therefore the decay rate 
is independent of /i and \x. 

We can compare our result with the differential decay rate for B — > X s ^ in the endpoint 
region at leading order [3| , 

dY G 2 m 4 

(B -> X s >y) = — b x z aH 1 {m b ^ Q )S{x 1 , // ), (40) 



dE^ " S,J 16vr 4 7 



where x 7 = 2E^/m b , and a is the fine structure constant. if 7 is the hard coefficient 

HJm b ,(io) = |V t6 V t :| 2 |C 7 | 2 |d + C 2 | 2 , (41) 



with C1 + C2 = l+0(a s ). Here we have used the operator basis suggested in Ref. [29j for the 
Wilson coefficients, which is equivalent to the one in Ref. [ill ]. Note that 5, the convolution 
of the jet function and the B meson shape function, appears exactly as in B — > XM. 
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Therefore if we take the ratio of these two decays, this factor cancels out, reducing the 
theoretical uncertainty. In the £77(3) limit, the ratio is given by 



dT(B^XM) dT{B^X sl ) 



dE 



M 



dE„ 



'J-'M- 



ami H^m^fio) 



(42) 



which is only a function of the hard coefficients (with Hj$ including the convolution with 
the LCDA). The ratio does not depend on detailed information about the 5-meson shape 
function. If charming penguins are present, this result is modified as discussed in the next 
section. 



V. CHARMING PENGUINS 

The size of the nonperturbative charming penguins in two-body B decays has been de- 
bated recently. Semi-inclusive decays B — ► XM can lead to new insight. Unlike two-body 
B decays, where additional nonperturbative parameters related to the B — > M form factors 
enter the predictions, the only nonperturbative parameters in Eq. (J42|) are the LCDA. If 
there are experimental deviations from Eq. (pf2"j) that exceed the uncertainties from sub- 
leading corrections when we compare processes with and without charming penguins (such 
as B — > X°(p), they would then unambiguously confirm the nonperturbative nature of the 
charm contribution. 

A typical charming penguin contribution is shown in Fig. 01 When the momentum 
transfer through the gluon is close to 4m^, the intermediate charm quark pair is almost 
on-shell and should be treated nonperturbatively, governed by usoft interactions. The long- 
distance contribution can be power counted as leading order in SCET [22] and cannot be 
disentangled from the bound state of the bottom quark. We can write the momentum of 
the charm quark pair as 2m c v^ c + /c M , where v^ c is the velocity of the charm quark pair, and 

is the residual momentum of order Aqcd- Note that v£ c is not the usual small velocity 




FIG. 3: A typical charming penguin contribution, with charm quarks in the loop, and a, (5 are 
the color indices. The outgoing momenta p and pm are n and n-collinear, respectively. Usoft 
interactions are not shown. 
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parameter v* in NRQCD. In the rest frame of the B meson, we can write 



2m c vg 




(43) 



where v 2 c = 1. The momentum fraction of u of the antiquark in meson M is given as 



where u = 1 — u and %m is close to 1 near the endpoint. 

There are two possible scenarios when we take appropriate limits of the charm quark 
mass compared to m^. First, we can take the heavy quark limit nib — > oo with m c ~ \JmbA. 
In this case, r 2 is of order A/ 'm& and the momentum of the charm pair, or the charm 
quark itself, becomes hard-collinear in the n direction because u = 4t 2 /xm is of order 
A/m&. Therefore the outgoing antiquark with the momentum upu is an usoft quark, and 
the exchanged gluon has offshellness of order y/m^K. The long-distance charming penguin 
contribution, shown in Fig. H|a), can be treated in the same way as that of the quarks 
with small mass, shown in Fig. b). Since the usoft-collinear interaction is suppressed at 
least by A ~ y/A/rrib, the leading long-distance contribution is suppressed by m c / nib at the 
operator level. In addition, the unbalanced endpoint configuration for the meson M gives 
the endpoint suppression of order A/rrib. Therefore, the long-distance contribution in this 
limit is suppressed by m c A/m 2 . This power counting is in agreement with the expectation 
that the contribution of the charming penguin in this limit gives the same contribution as 
massless quarks, which is suppressed by A 2 /ml |3C|. The reason why it is not of order 
rriqA/m 2 is because the quark mass insertion is replaced by the insertion of . 

The second scenario is to take the heavy quark limit nib, rn c — > oo with fixed r, which 
is motivated by the fact that u = At 2 /xm is near the central point experimentally. In 
this case, the power counting is different. The charm quark is regarded as a heavy quark 
and there is no endpoint suppression. The exchanged gluon in Fig. |3] has large offshellness 
Am 2 ~ ml, and is integrated out to obtain an operator of the form [cc] nr. [£n£n]- Here we 
suppress the Dirac structure and color indices, and the charm quark is treated as a heavy 



FIG. 4: a) Nonperturbative penguins where charm quarks in the loop are treated as hard-collinear, 
with an insertion of C mc of order A in SCETj. b) one of the nonperturbative penguin contributions 
for massless quarks with the insertion of and , each of which is suppressed by A in SCETj. 
Dashed line is an usoft quark, leading to additional endpoint suppression in forming a meson M. 



4r 2 = ux M , 



r = m c /m b , 



(44) 
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quark described by NRQCD. The nonperturbative charming penguin contribution is then 
obtained from the time-ordered product of this operator with the four-quark operators in 
the weak Hamiltonian of the form [£ n ^] [cc]nr- This leads to a contribution ~ a s {2m c )v*T C c 
in agreement with Refs. 19, 2^|, where the nonperturbative function T C c is comparable in 



size to the leading-order shape function /. Furthermore there is no endpoint suppression 
because the outgoing antiquark forming the meson M is collinear with momentum fraction 
u = 4t 2 /xm, corresponding numerically to the the central region in the LCDA of the meson 
M. 

We think that the second scenario is more plausible based on the actual value of m c /m h . 
As this is the more conservative of the two, we consider the nonperturbative charming 
penguins in the second scenario, which give larger contributions than the first scenario. As 
explained above, in the heavy quark limit m^ c oo with r fixed, the charming penguin 
can be of leading order, which can be expanded in powers of A/m^c and a s (2m c ) in a 
consistent way. In this scheme, at leading order in l/m;,. c , the contribution is factorized into 
the n-collinear part, the jet function, and a nonperturbative function. The derivation of the 
factorized form is presented in detail in Appendix El At leading order in A/m& iC , and to 
first order in a s (2m c ), the contribution of the charming penguin to the differential decay 
rate can be written as [see Eq. (|B20|) 

dT 5c (B->X q M) = G| /Mx 3 f A (,) as(2mc) c bm u ( u = 1 _^l\. 2 ReT^, (45) 

where c q M is defined in Appendix |B] and T^f is the hard kernel given in Table |I] Here 
M = P, Vl, while the contribution of Vr is l/mfe c suppressed because of the spin flip. The 
function T C c does not depend on the outgoing meson M or the jet X because is given 
by the nonperturbative effects arising only from the usoft interactions of the on-shell charm 
quark pair in the B meson. Hence, up to the B meson flavor, T C f is universal in all the 
decay modes where charm penguins contribute, and its size is experimentally measurable 
from various decay modes. In the isospin limit the T C c functions in B and B~ decays are 
the same, and are equal to T C c in B s decays in the SU[3) limit. Due to its nonperturbative 
nature it can, however, have a nonzero strong phase |l9(. 

In summary, the differential decay rate for the semi-inclusive decays at leading order in 
l/rrib, including the nonperturbative charming penguin at LO in l/m c and a s (2m c ), is 



dV , „ „ G 



B - X q M) = -f / 



; 2 mix 3 . 



dE M o7T 



H$(m b ,iio)S(xM,tM}) ( 46 ) 



+a s {2m c )\^c B q M <p M (l - —) 2ReT^*T 

V X M J 

Phenomenological implications of this expression will be discussed in section I VI II 
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VI. ESTIMATES OF SUBLEADING CORRECTIONS 



To predict the branching ratios for B — ► XM more accurately, Eq. (J46|) can be systemat- 
ically extended to higher orders in l/m& and a s . In this way one could also unambiguously 
determine whether a possible future discrepancy between experiment and predictions based 
on Eq. (J42)) is due to nonperturbative charming penguins or higher-order corrections. A full 
analysis of the higher-order corrections is beyond the scope of this paper, but we identify typ- 
ical subleading corrections and estimate their size. The subleading corrections, suppressed 
by powers of A/m&, are of two types: (i) corrections to the heavy-to-light current leading 
to the subleading S-meson shape functions, some of which are already well known from the 



351 ]. and (ii) corrections 



analyses of B — > X s ^ and B — > X U W inclusive decays jjl. 32. 33. 34. 
to the n-collinear currents forming the light meson M, which appear as the twist-3 LCDA 
and the 577(3) breaking effects in the twist-2 LCDA. 

Let us consider the corrections of the first type. The convolution S of the jet function 
and the B meson shape function in Eq. (|37|) can be expanded to higher orders in l/m& as 



S{x M ) = S®{x M ) + (sfi>(x M ) + Sg\x M j) + ■ 



(47) 



where «s£j is the subleading corrections to the heavy-to-light current and S^ 1 is the usoft 



(i) 



corrections to the n-collinear current. Using the results of Ref. [35| , the sum can 
be related to the imaginary part of the time-ordered product, W^, 

1 

m b 

where the factor n tl /2 comes from the n-collinear current (£n7 M £n) = ^ M (6r^£n)/2- Taking 
the ratio with respect to B — > Xg'y gives the difference 



■ (5W) + 5 g)) = ww ! 



(48) 



dn ■ I v(n ■ I) 5(mf,(l — 



(49) 



+ n-l) = Ml 



x 



nit 



where x 



is chosen. The subleading shape functions v and H 2 are de- 



fined in Refs. 35] and |31j], respectively. In particular, H 2 is proportional to A2 = 
(B v \b v g s a ^G^ b v \B v ) / 12 ~ 0.12 GeV 2 . Taking a broad cut E M > 2.0 GeV, this contri- 
bution should not exceed 10% compared to the leading contribution, unless there is an 
enhancement in the coefficient. 

The subleading correction <S^ comes from the usoft interactions with the n-collinear 
currents, which lead to new subleading B meson shape functions. The subleading operators 
obtained by inserting the usoft gauge-invariant term Y^ip^ s Yn are suppressed by A 2 in 
SCETi, but they should be included in SCET n because they are suppressed by A/m&. The 
nonzero contributions come only from the color-octet four-quark operators 

2(q n Y^P L T a YX) ■ [q'^P L , R T a Yii]% s Y w -^ji q £ 

1 



{ -Ku) 



2(q n Y^YP L T a YX) ■ q^i—^Y^p^T^P^q^ 



(50) 
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where the flavor and chirality structure is the same as the corresponding leading operators 
Oi in Eq. (jSJ). Because of reparameterization invariance jH^, the Wilson coefficients of 
Of\(u) + Of B (u) are the same as those for the leading color-octet operators in Eq. ()12|). 
which are presented in Appendix 1X1 

After some calculation, the matrix elements of these operators which do not vanish triv- 
ially from the flavor content are nonzero only for specific values of i due to the helicity 
structure. They are given by 



{Ofl{u)) 



ifM <Pm(u) 



N 



{OfUu)) = (5 



VM 



U 



S 



(X q \q n Y^ S Y W P l YX\B), (^ = 1,2,4,5), 



(51) 



PM j 



ifjM 4>m(u) 



N 



u 



(X g \q n Y^YJYitpi]p L b v \B), (i = 3,6), 



where 5pm, <Vm are Kronecker deltas, and we use Eqs. (|B10Jl and (|B11J) in Appendix El for 
the matrix elements of the collinear current. For the matrix elements of the heavy-to-light 
current, we need the time-ordered products with J H , 



T 



(2) 



2i 



A 



T, 



(2) 



m h m B 
2i 



d'z e-^(B\Tf H (z) q n Y^Jft^ P L Y^b v {Q)\B) , 



(52) 



B 



d 4 z e-^ z (B\Tf H (z) q n YiY«\Yiipi 



B). 



m b m B 

They can be factorized into the jet function and subleading 5-meson shape functions, 



T 



(2) 



A.B 



-m h (l-x M ) 

where the subleading shape functions are defined as 



dn ■ I f^ B (n ■ I) Jp (m b (l — i M ) + n ■ I + iej . 



f [ I\n-l) = — (B v \b v YJ{n-l-n-id)'^ Y^Y n Yib v \B v ), 
rrib 



T 



(53) 



(54) 



n-l) = —(B v \b v Y n 5{n-l-n-id)'^Y^ 

e(l) 



Y^i . 



b v (0)\B v ). 



The shape functions f A B are different from the subleading shape functions appearing in B — > 
X s 7 and B — > X U W, due to the presence of Y^\ which cannot be neglected at subleading 
order. At present we cannot estimate their size, but there is no reason to expect a dramatic 
enhancement from the insertion of Y^\ However, these contributions can be numerically 
significant in the decays that are very small at LO in These are the color-suppressed 

AS = tree decays B s — > {n , p°, the QCD penguin-dominated AS = and AS = 1 

decays B — ► (f)X°, B~ — ► <\>X~ and B s — > loX^ and the A« ' part of the amplitudes in 



B° s — > {rr°, p°}X® s . For these decays, the Wilson coefficients of the LO operators convoluted 



with the asymptotic LCDA (\(f> ® C$\ ~ 0.08 and \<j> ® (C 5 + C 6 

1.9 and I 



•> 0.005) are much smaller 

than those for the color-octet operators (\(f> (g> C%\ ~ 1.9 and \<j) <g> (C 5 + C 6 )| ~ 0.15). The 
subleading contributions can thus be numerically large in spite of the 1/2N suppression. 
Note that this is not a sign of failure of the l/m& expansion, but due to the hierarchy of the 
Wilson coefficients. 
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FIG. 5: Subleading usoft interactions induced from the n-collinear currents. In Diagrams b) and 
c), the dots denote the A-suppressed interaction in . 

There is another contribution shown in Fig. from the time-ordered products of the 
n-collinear currents and , 

4 } = hWn(Y4il[kY v ) J^wkflitfr + y^t Wv J- (YiifaY*) tufa. (55) 

The intermediate legs in Fig. are hard-collinear and give a jet function of the form 1/n-k. 
The relevant LCDA for M are suppressed by A/m b due to the presence of iD±. But it is 
not known whether this process factorizes, and we leave a full analysis for future work. 

We now consider the contributions from the (S — P) x (S + P) four-quark operator. At 
leading order it matches onto q' n (l — 75)6 • %(1 + 75)?^, which vanishes because of spin 
symmetry. At subleading order it matches onto 



Of +P = -2 J2 (q'ntPRY r %)(q w P R ^W^W JT 4q' fl ) + 

q'=u,d,s 



(56) 



where for the heavy-to-light current we have used the relation (2v M = n fl + n M ) 

2i n P L YX = 2q' n P L ^YX = q' n fP R YX- (57) 
In semi-inclusive B decays, the amplitude from this operator factorizes using the twist-3 



LCDA 38J. At order l/m b it contributes through the time-ordered product with the leading 
heavy-to-light current as 

T s+P = — I ' d A z e -*"-*<fl|T4(z) q'jiP R YX(0)\B}. (58) 

m B J 

This vanishes because 

(B\Tf H (z) qjtP R YX(0)\B) oc (B\b v Y n (^)fP L ifP R YX(0)\B) = 0. (59) 
Therefore the nonzero contribution comes from the time-ordered products of Og} p with the 



subleading operators of suppressed by m q /m b 



29|, or from the time-ordered products 
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of { p +P with itself. Both are of order l/mf, but the latter may not be numerically negligi- 



ble. In the QCD factorization approach |2J|, the contributions corresponding to C^+p lead 
to formally suppressed but numerically large "chirally enhanced" contributions. In SCET 
®S+p * s a ^ so f° rma lly suppressed, while its matrix elements remain unknown and could be 
numerically large. Because of this uncertainty, the decay rates and the CP asymmetries 
presented in the next section for modes with small tree-level amplitudes should be regarded 
only as a rough estimate. 

Finally, we discuss the SU(3) breaking corrections to Eq. The SU(3) breaking due 

to different light-quark flavors in the inclusive jet is suppressed by rn^jm^h. [2^ and thus 
negligible, but the 577(3) corrections due to the strangeness content of meson M are not 
negligible. These are realized in SCET by inserting the strange quark mass term 



n ■ iD~ 



2^' 



(60) 



in the leading n-collinear currents with the final-state s quark in Fig. |3 It can be written 



as 



AJu) 



■(M\ T [s n 4P L q n (0)] u -t / d 4 zC m (z) |0). 



(61) 



I'm Em 

For the final-state s quark, a hermitian conjugate of Eq. (JBTj). A s (u), is needed with A s {u) = 
A 8 (u). If m s is comparable to iD± ~ A, C m is of leading order in SCET n , suppressed only 
by m s /A. 

The SU(3) breaking affects meson decay constants and the LCDAs. One finds to leading 
order in S , f/(3)-breaking <Pk(u) = (f>(u) + A s (u), 4>k(u) = 4>(u) + A s {u) and 4> v (u) = 4>{u) + 
2(A s (u) + A s (u)) /3 for the LCDA of K~(K°), K + (K°) and rj, where <f>{u) is the pion LCDA. 
Since 4>M(u)du = 1 = L 4>{u)du one has a constraint 



A 8 {u)du = 0. 



(62) 



It is also straightforward to check that these LCDA satisfy the relation 41] 4> n (u, /i) + 



3^(w,/x)=2 (j) K +(u,n) +(j) K -(u,n) =2 4> K o(u,jj) + (j) R0 (u,n) 



From power counting one 



expects the relative size of the SU(3) breaking contribution to be of order m s /A ~ 20% 
Recent QCD sum rules predictions can be found in Refs. 



42 



44, 





FIG. 6: The insertion of the strange quark mass in the n-collinear currents in SCETn- Diagrams 
a) and b) represent the possible strange quark mass insertions for the strange quark in Eq. H61|) 
where the external particle is a strange quark. 
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VII. PHENOMENOLOGY 



In this section we collect the predictions for B — ► XM decay rates and direct CP asym- 
metries, defined as 

A = dT(E -> XM)/dE M - dT{B -+ XM)/dE M 
CP dT(B -> XM)/dE M + dT(B -> XM)/dE M ' 

while treating the charming penguins as perturbative. Once the experimental data become 
available, one of the modes can be used to determine the nonperturbative charm contribu- 
tion Tec and then modify the predictions according to Eq. (|46)l . To reduce the hadronic 



uncertainty, we normalize the branching ratios of 5" — > XM, B — > XM, and B° s — > XM 
to the decay rates — > X°7, — ► X~7, and £> s — > X° 5 7 in the endpoint region, re- 
spectively [See Eq. f!42|) .] The only remaining nonperturbative input is then the light meson 
LCDAs. Expanding in terms of the Gegenbauer polynomials, 

oo 

Mu,fi) = &uu[l + Ys<(v)C 3 n /2 (2u - 1)], (64) 

71=1 

we truncate the series at n = 2 and use isospin symmetry to set af 1 = for mesons not 
containing a strange quark. We fix the remaining coefficients using the results from QCD 
sum rules, while conservatively doubling the errors quoted in the literature. This gives at 



H = 2 GeV: af = 0.05 ±0.05, af = 0.23 ±0.23 [43], < = 0.08 ±0.13 



441, al = 0.09±0.15 



441, while for the u LCD A 



45J, a£ = 0.07 ± 0.08, a p 2 = 0.14 ± 0.15, a v 2 = 0. ± 0.15 
a 2 = 0. ± 0.2 is used for lack of better information. 

Direct CP asymmetries in Eq. ()63|) are nonzero only in the presence of nonzero strong 
phases. These are generated by integrating out on-shell light quarks in a loop when matching 
full QCD to SCETi at NLO in a s . We therefore use the NLO matching expressions for the 
Wilson coefficients Cf at \i — even though the evolution to the hard-collinear scale 
/io ~ V 'Awifc is performed at NLL. Note that this running cancels to a large extent in the 
ratios of decay rates (only the running of a„ (//), n > 1 remains), giving in effect the Wilson 
coefficients with NLO accuracy at the hard-collinear scale Hq. 

For definitiveness, we choose /io = 2 GeV for the hard-collinear scale, which corresponds to 
the experimental cut p 2 x < (2 GeV) 2 on the inclusive jet invariant mass. The corresponding 
predictions are listed in Tables HP and ITTT1 for AS* = 1 and AS = decays respectively. The 
predicted partial decay widths dT(B — > XM) j dEu in principle depend on the light meson 
energy Em- In the endpoint region the dependence on xm = 2Em/™>b = l±( m M ~ Px)/ m2 B 
is, however, a subleading effect, 2 so we set %m = 1 in Tables ITT1 and IIII1 

The two errors quoted in Tables |n] and Mil are an estimate of subleading corrections and 
the errors due to coefficients of the Gegenbauer expansion of LCDAs. Since the predictions 



2 Numerically, for p 2 x < (2 GeV) 2 one has E w > 2.26 GeV compared to m B o/2 = 2.64 GeV. The same 
p 2 x cut corresponds to higher Em cut for heavier mesons, for instance, for (j> mesons the same cut on p x 
corresponds to > 2.36 GeV. Thus ibb/2 — Em ~ A with 1 — xm ~ 0(A/ms). 
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are made to NLO in a s (rrib) but only to LO in l/m&, the largest corrections are expected to 
arise from the l/m& terms. These are estimated by independently varying the magnitudes 
of the leading terms proportional to A^'j^ by 20% ~ 0(A/m b ) and the strong phase by 
5°. This latter variation estimates the error on the strong phase arising from the uncalcu- 
lated a s {mb) /rrih or a 2 (mt) terms. A 100% error is assigned to predictions for branching 
ratios in color-suppressed tree and QCD penguin-dominated AS = decays where the l/m^ 
corrections are sizable compared to the leading results due to the hierarchy of Wilson coef- 
ficients. No prediction on CP asymmetries is given for these modes or for the affected QCD 
penguin-dominated AS* = 1 decays. 

To understand better the relative sizes of different branching ratios, it is useful to split 
the amplitudes for the semi-inclusive decays according to the CKM elements. Using the 
unitarity of the CKM matrix A^ = — Xu^ — Ac , the amplitude can be rewritten in terms 
of the "tree" amplitude T b ^mx and the "penguin" amplitude Pb-^mx as 

A(B - MX) = X^T B ^ MX + X^Pb^mx, (65) 

with q = d,s for AS = 0, 1 decays respectively. The "tree" amplitudes receive contributions 
from m Eq. (JBJ), the "penguin" amplitudes from 0\ 2 (charming penguins), while the 
QCD and electroweak penguin operators contribute to both amplitudes. The combinations 



Mode 




Br(Mode)/Br(S X s j) 


Exp. (2-body) 


A C p 


B°^ 


K~X+ 


0.16 ±0.09 ±0.05 


> 0.078 


0.30 ±0.16 ±0.01 


5°-> 


K*-X + 


0.28 ±0.16 ±0.07 


> 0.026 


0.31 ±0.16 ±0.02 


U° — > 




0.22 ±0.13 ±0.03 


> 0.034 


0.0089 ± 0.0050 ± 0.0016 


Br 


■K°X- 


0.20 ±0.11 ±0.06 


> 0.067 


0.0097 ± 0.0048 ± 0.0006 


B- 


■K*°X- 


0.34 ±0.19 ±0.08 


> 0.040 


0.0084 ± 0.0046 ± 0.0019 


B- — l > 


■ <px- 


0.22 ±0.13 ±0.03 


> 0.035 


0.0089 ± 0.0050 ± 0.0016 






(1.0 ±0.6 ±0.2) x 10~ 2 






BS- 


K'Xf 


0.16 ±0.09 ±0.05 




0.30 ±0.16 ±0.01 






(2.4 ± 1.4 ±0.5) x 10~ 2 








uX% 


(2.8 ±3.4 ±0.7) x 10~ 3 






«N 


K*-Xf 


0.28 ±0.16 ±0.07 




0.32 ±0.16 ±0.02 



TABLE II: Predictions for decay rates and direct CP asymmetries for AS = 1 semi-inclusive 
hadronic decays are given in the second and fourth column, respectively. The first errors are an 
estimate of the l/m& corrections, while the second errors are due to errors on the Gegenbauer 
coefficients in the expansion of the LCDA. The third column gives lower bounds on inclusive decay 
rates obtained by summing over measured two-body decays and normalizing to b — ► 57 decay with 
E min = 2.0 GeV (90% CL lower bounds are used). 
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Mode 




Br(ModeJ/lJr(i3 — ► A s 7) 


bxp. (2-body) 


A C p 


— 

B 


7T-X+ 


0.67 ±0.37 ±0.14 


> 0.038 


-0.040 ± 0.021 ± 0.004 


— n 

B -»■ 




(9.1 ±5.3 ±3.1) x 10~ 3 


> 2.0 x 10~ 3 


-0.15 ±0.11 ±0.01 


s -► 




(2.0 ±2.0 ±0.1) x 10~ 4 






-r=:0 
5 -► 


p -x+ 


1.76 ±0.97 ±0.35 


> 0.10 


-0.039 ± 0.021 ± 0.004 


— 

B -> 


K*°X® 


(1.4 ±0.8 ±0.5) x 10~ 2 




-0.17 ±0.11 ±0.03 


B~ — 5 


■ K°X- 


(9.1 ±5.3 ±3.1) x 10~ 3 


> 2.5 x 10~ 3 


-0.15 ±0.11 ±0.01 




■ cf)X- 


(2.0 ±2.0 ±0.1) x 10~ 4 






5" -i 


■ K*°X~ 


(1.4 ±0.8 ±0.5) x 10~ 2 




-0.17 ±0.11 ±0.03 


— 




(4.1 ±4.1 ±2.6) x 10~ 3 




— 


3- 


7T-X+ 


0.67 ±0.37 ±0.14 




-0.040 ± 0.021 ± 0.004 






(1.3 ± 1.3 ±0.7) x 10~ 2 








p-xf 


1.76 ±0.97 ±0.35 




-0.039 ±0.021 ±0.004 




uXl 


(1.1 ±1.1 ±0.9) x 10~ 2 







TABLE III: Predictions for decay widths and direct CP asymmetries of AS = semi-inclusive 
hadronic decays. The first errors are an estimate of the l/m^ corrections, while the second errors 
are due to errors on the Gegenbauer coefficients in the expansion of the LCDAs. The third column 
gives lower bounds on inclusive decay rates obtained by summing over measured two-body decays 
and normalizing to b — > 57 decay with E m \ a = 2.0 GeV (90% CL lower bounds are used). 

of the CKM elements exhibit the following hierarchy 

A^-A 2 ,, -A 3 ,, A^~A 4 C , (66) 

where Ac = 0.23 is the Cabibbo angle. In AS = decays, the two CKM factors in Eq. (}65|) 
are of comparable size. In AS = 1 decays, on the other hand, there is a hierarchy between 
the two terms in Eq. (J55J) since |Ai /Ac | ~ \ 2 C . To first order in this small ratio, the 
quantity 

A^iB^ MX) = -2I m (|) I, (^) , (67) 

with — 2Im(Au ^/Xc ) = 0.037, which sets a typical size of the CP asymmetries. The size 
of the direct CP asymmetries also crucially depends on the ratio of "tree" over "penguin" 
amplitudes, as can be seen in Table [H] This can be estimated from the sizes of the Wilson 
coefficients at fi — mi, (convoluted with the asymptotic form of LCDA) that are given in 
Table IIVI The modes that receive contributions from the operator Of, B s — > Xf 
and B KW-X+, have T B -*mx > Pb^mx and thus have larger CP asymmetries. The 
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C i 


pu 
u 2 


pc pc p 
L 1 L 2 C 3 


pu 


nc 
u 4 


C5 Cq 


Abs 


0.89 


0.080 


0.0011 0.012 0.00037 


0.029 


0.032 


0.010 0.0061 


Arg 


0.9° 


-99° 


79° 181° 7.5° 


-150° 


-163° 


14° -151° 



TABLE IV: The magnitudes and strong phases of the Wilson coefficients at fx = mj, (using the 
notation £3,5,6 = C3 5 6 = ^3 5 6) convoluted with the asymptotic form of the LCDA cj) = 6uu. 

rest of the modes listed in Table HU do not receive these large tree contributions and thus 
have smaller CP asymmetries. Note that the direct CP asymmetries are nonzero only if the 
two interfering amplitudes in Eq. (f63j) have different strong phases. In the AS = decays 
B° — > (f>X° and B~ — > (pX~~ , the two amplitudes Tb^mx, Pb^mx are the same at LO in 
1 / mb and the CP asymmetries vanish. This may change at higher orders, but no prediction 
for Aqp for these modes is given in Table IIII1 

For color-suppressed two-body decays, the leading order tree amplitude in SCET : comes 
from O2. However, when matching onto SCETn, the hard-spectator contribution from 0\ 
can compete with the leading order term. For semi-inclusive decays considered in this 
paper, in which the spectator quark does not enter the outgoing meson M, there are no 
hard-spectator interactions. Thus, due to the hierarchy of the Wilson coefficients (using 
values in Table UVl 



f du6uuC%{u) -0.10 f duQuuCKu) , (68) 
Jo Jo 

numerically (lH l2o| . semi-inclusive tree amplitudes receiving contributions from O2 are 
smaller than the tree amplitudes due to 0\. The color-suppressed tree decays are then 
more sensitive to 1/rrib corrections, as discussed in the previous section. These may be 
especially important for the decay B° s — > wX s ° in which a cancellation between different 
contributions occurs for central values of input parameters. A strong dependence of the 
predictions on a 2uJ is thus found with Br(B° s uX?)/Br(B° s 7 X S ) G [0.003,0.027] 
for a 2aJ G [—0.3,0.3]. Sizable 1/mb corrections are expected in all the modes without the 
charming penguin contributions: B —* <f)X°, B~ —* <pX~ , and B s — > MX S ° 5 (M = 7r ,p ,a;). 
These decay modes are a good experimental source to analyze the corrections at order 
1/mb. A testing ground for the charming penguins are the processes in which the tree- 
level amplitudes are not suppressed, and there is a charming penguin. They correspond to 
processes with C\ and Cf in Table |U 

In Tables |H] and IIHI we also give the experimental lower bounds on the predicted semi- 
inclusive branching ratios. These were obtained by summing over the already measured two- 
body decays and normalizing them to Br(6 — > S7) = (317±23) x 10~ 6 for E 1 > 2.0 GeV. The 
two-body channels for which only upper bounds are known were not used in the estimate, 
nor were the decays to more than two hadrons in the final state. 

Experimentally, the semi-inclusive hadronic decays can be measured either by summing 
over exclusive decays or by performing a truly inclusive measurement where only the flavor 
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and charge of the decaying B meson and of the isolated energetic light meson M are tagged. 
For these measurements a first step might be made by making an even more inclusive 
measurement where only the flavor, but not the charge of the initial B meson is tagged. 
Theoretically simple predictions can be made for B~ fW -> K S:L X S , B~ fW -> K*°X S 
and B~ /B~° — > <pX s decays, where B~ /B~° denotes a sum over the decay widths, T(B~ — > 
MX) + T(B° — > MX). Using isospin symmetry, the following relations hold in the endpoint 
region due to factorization at leading order in l/ra&: 

Br(B- - K S , L X;) = Br(£° - K S , L X° S ), Br(B~ - K*°X;) = Br(B° - K*°X° S ), 
Br(B- -> <j>X~) = Br(B° -> 0X S °), Br(S~ -> = Br(5° -> 0X°), (69) 

so that 

Br(B-/B° - K 5 , L X S ) = l -Bi{B-/B^ - K°X S ) = ^Br(£T - K°X7), 
Br (5-/5° -> iT°X s ) = ^Br( J B- -> iT°X;), (70) 
Br(5-/5° - 0X (S) ) = l -Br{B~ -> 0X ( " } ). 
For the direct CP asymmetries of these more inclusive modes, we find 

A CP (B-/B° - K S , L X S ) = A CP (B-/B^ - K°X S ) = A CP (£T - K°X;), 
A CP (S-/S° - 7T°X S ) = A CP (iT - 1T A7), (71) 
A CP (B-/B^ - 0X W ) = A CP (£T - 

and A CP (B-/B^ -> 0X u+d+s ) ~ A C p(£T/£° -> 0X S ). 

Furthermore, for S — > 0X and B K X decays an even more inclusive measurement 
can be made, where the strangeness content of the inclusive jet need not be determined, sim- 
plifying the measurement. Since Br(£>~ / B — > 0X S ) ^> Br(B~ /B — > 0X), the theoretical 
prediction for this inclusive measurement is Br(B~ /~B° — > 0X u+d+s ) ~ Br(£?~/_B° — > 0X S ) 
valid up to corrections at the percent level. A similar simplification occurs in B~~ / B — > 
K*°X S decays, since the decays B~ — > fT*°X~ and 5° — > fT*°X are absent. Therefore the 
strangeness of the inclusive jet is fixed automatically and need not be determined experi- 
mentally. An important part of the measurement is that the flavor of K*° is tagged from 
the decay K*° — > K + n~ . On the other hand in B~ / B — > Ks,lX s decays, since there are 



contributions with the spectator quark ending up in K" from B" — > i^°X°, the strangeness 
content of the inclusive jet should be determined from experiment. 



VIII. CONCLUSIONS 

In the framework of SCET we have considered semi-inclusive, hadronic decays B — > XM 
in the endpoint region, where the light meson M and the inclusive jet X with p 2 x ~ Am?, 
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are emitted back-to back. We have considered the decays in which the spectator quark does 
not enter into the meson M. In SCET the four-quark operators factorize, which allows for 
a systematic theoretical treatment. After matching the effective weak Hamiltonian in full 
QCD onto SCETj, the weak interaction four-quark operators factor into the heavy-to-light 
current and the n-collinear current. The forward scattering amplitude of the heavy-to-light 
currents leads to a convolution S of the jet function with the I?-meson shape function, while 
the matrix element of n-collinear currents gives the LCDA for the meson M, leading to a 
factorized form for the decay rates. The two nonperturbative functions, the convolution S 
and the LCDA, are the only nonperturbative input in the predictions for B — > XM decay 
rates at leading order in 1/m.ft. Furthermore, the same convolution S appears in B — > A s 7 
decay and drops out in the ratio of B — > XM to the B — > X s 7 rate and in the prediction 
for direct CP asymmetries. 

This greatly reduces hadronic uncertainties, since the remaining nonperturbative input, 
the LCDA, is well described by its asymptotic form, corrections to which can be obtained 
from other experiments or from QCD sum rules. The Wilson coefficients can be perturba- 
tively computed and are then evolved to the scale /x ~ y/Am b using the NLL expressions. 
In the ratios the multiplicative RG evolution factors almost cancel. The predictions for 
branching ratios and CP asymmetries are then given at NLO in a^mj) and at LO in 1/m.f, 
and are collected in Tables ITT1 and IIII1 Numerical values are given in the limit of perturbative 
charming penguin due to a lack of experimental data, while the formalism used is extended 
to the case of nonperturbative charming penguins. To leading order in SCET, the charming 
penguin contribution factorizes and is given by a universal nonperturbative function jF cS 
describing the usoft interactions between the on-shell charm pair and the bound state of the 
b quark. In particular, does not depend on the final meson M or the flavor content of 
the inclusive jet, but only on the flavor of initial B meson. 

We have also estimated subleading corrections and identify potentially large subleading 
usoft contributions coming from the n-collinear sector giving rise to color-octet operators. 
These contributions can be of appreciable size, compared with the leading contributions 
when the leading contributions are suppressed by Wilson coefficients. This, for instance, 
happens for color-suppressed tree decays and QCD penguin (not charming penguin) dom- 
inated decays. Other contributions such as the (S — P) ® (S + P) operators that have 
been argued to be large in exclusive B decays, on the other hand, vanish to first order in 
l/m^, but are present at higher orders. Similarly, subleading corrections coming from the 
heavy-to-light sector and giving subleading 5-meson shape functions largely cancel in the 
ratio with the rate B — > X s j. 

In conclusion, semi-inclusive hadronic B decays are a good test field to clarify many 
hadronic uncertainties common to two-body exclusive B decays and the inclusive B decays 
at the endpoint. The factorized results provide us with a simplified view on the diverse chan- 
nels of hadronic B decays and enable us to consider them rigorously within the framework 
of SCET. By investigating decays without charming penguins, we can test whether the for- 
malism is working. Then by looking at modes where the charming penguin can contribute, 
we can potentially see whether or not the charming penguin give a large contribution to the 
decays. 
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APPENDIX A: THE WILSON COEFFICIENTS AT NLO 



The matching of the weak Hamiltonian in Eq. flj|) from full QCD to SCETi was calculated 
at NLO in a^m^) first in Refs. [21|, and then in Ref. For the detailed matching 

procedure in obtaining the Wilson coefficients, the reader is referred to Ref. 0|. Here we 
translate the results to the basis choice of Eqs. (jHJ) and (|T2*j) . The Wilson coefficients for 
operators 
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CUv) = 8, 



up 
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with the shorthand notation JC(v) = —6 — 7T 2 /12 and 
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W An 
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7T 

12 
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\nv 



- 2 In 2 v - 2Li 2 (l -v)-(v<->v 

~lnv \nv~ 



F{v) =F(v) + Qin + 24 + 3(v - v) 

IV V 

The contribution of a fermion loop and the gluonic operator to C\{v ) is given as 
L p =|(c7 a + 2C 3 + 5C 4 - C 9 + - (C 3 + 3C 4 + 3C 6 - y)(?(0) 

- (c 4 + C 6 + C 8 + C w )g(z c ) - (C 3 + C 4 + C 6 - ^(C 8 + C 9 + Cio))G(l) 

- |(c 5 + C g - l -C 7 ) - C 1 (5 up G(0) + 8 cp G(z c )), 



(Al) 
(A2) 

(A3) 
(A4) 



(A5) 
(A6) 



(A7) 



Note that C^ 5 6 (v) — Cf 5 6 (u), so we also use the notation £3,5,6^) = Cf 5 6 (v). 
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where Zf — rnj/ml and 



—4 / dw w(l — w) In Zf — w(l — w)v — it 



G(z f , 



The Wilson coefficients for the octet operators in Eq. ()12j) are 
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APPENDIX B: NONPERTURBATIVE CHARMING PENGUIN IN THE 
HEAVY QUARK LIMIT 

In this Appendix we show that in the heavy quark limit, m;,,m c — > oo with r = m c /mb 
fixed, the charming penguin contributions to the decay rates factorize in SCET into hard, 
jet, collinear, and soft parts at LO in l/m C) 6. A typical charming penguin contribution is 
shown in Fig- El When the momentum transfer in the gluon is close to 4m 2 , the intermediate 
charm quark pair is nearly on-shell and can have usoft interactions. In the B meson rest 
frame, the velocity of the b quark can be written as v M = (n^+W) /2 with = (1, 0, 0, 1) and 

= (1, 0, 0, —1). In this frame, the on-shell charm quark pair has momentum 2m c v^ c + 
where k 11 ~ A is the residual momentum, while v£ c is the velocity of the charm quark pair 
with v 2 C = 1. It is given by 



1 n" 
2rY + 



uxm n 11 



2r 2 

where 4r 2 = uxm, with xm close to 1 and u—l — u 



1 



(Bl) 
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The charm quark pair annihilates into a gluon with off-shellness of order 4m\ 



mi. 



Integrating out the intermediate off-shell gluon gives a four-quark operator at leading order 
in l/m c 

O c5nn = ^n,.rT a qn^) (c. v ^,T a c v J, (B2) 
q 

where the charm quarks are treated as heavy. The collinear quark fields q n and q n are defined 



as 



5{u x -n-V)W^ 



(B3) 



where W n (Wn) is the collinear Wilson line in the n (n) direction from integrating out off- 
shell heavy charm quarks. Note that these collinear Wilson lines are the same as those 
from the h eavy b quarks in Eq. (fTTj) . This is a manifestation of Type-Ill reparameterization 
mvanance which states that the SCET Lagrangian and the collinear Wilson lines 

are invariant under — > n^/ f3 and n a — > f3n a with /3 close to 1. For example, the collinear 
Wilson line W„ is invariant under this transformation as 



n'-A 



w-v 



pn-A 
f3n-V 



n ■ A 
n-V 



(B4) 



which also holds for W n . This corresponds to the Lorentz invariance under a boost with 
(3 = mb/(2m c ) in the z direction, corresponding to transforming to the B meson rest frame. 
The usoft interactions can be decoupled from collinear interactions by introducing the usoft 
Wilson lines Y n and Y n and redefining the collinear fields [13]. This gives 



O c - cnn = ^n^T^) (c^Y^T^c^ 



(B5) 



where the collinear fields from now on will denote the redefined fields. The operator O c5nn 
satisfies gauge invariance in SCET 13, |46J, and the sub leading corrections to this operator 
will be of order A/m c ^- 

Let us next discuss the form of the nonperturbative charming penguin contributions that 
arise from the time ordered product of O c5nn with the operators in the weak Hamiltonian. 
We work out the details for the operator 0\ = 4{q^ v PLc){cr) u PLb) (q = d, s) that matches 
onto the SCET operator 



qbcc[^2) 



(B6) 



where a, j3 are color indices. The treatment of other operators is similar. The matrix 
element for the contribution of weak operator 0\ in Fig. El is then 



{MX\H§\B) = i f duduid^f d A y ^ . 

' ' Pn,Pn 



i{—2m c Vcc+Pn—Pn)-y 



(B7) 



(u,ux,y) O qbc5 {uj2 , 0) I B) , 
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where H C c n n = —4ira s / (uloi) and H q b C c = 1 + 0(a s ) are the Wilson coefficients of the 
operators C mn and O q b C c in Eqs. (|B5J) and (|B6j) . 

Using the factorization of n-collinear quarks from usoft and n-collinear degrees of freedom 
the matrix element (IB 71) can be rewritten as 



G/? \®c% M h M [ du 5(u-—)<f> M (u)H(u,x M ) (X\Q c5 \B), (B8) 

JO V X M ' 



V2N 



where 



(X\Q C5 \B) = -Ai J d 4 y(X\Tq n (y)rr M rPLT a (Yic. v J(0) 

x (c^ c Y^T a Y^ c ) (y) {c VB Y nlv P L YXm I B) . (B9) 



The delta function in (|B8|) is obtained from the exponent of the label momenta in Eq. (JB] 
using n ■ p — m b , n ■ p M = xu^b- The hard kernel is H(u, xm) = 47ia s (2m c )Ci/ (uxmitiI) + 
0(a1). In obtaining (|B8J) the relation 



(M\(q w ) a [q w 8(u-n- V^/n • Pa/)] JO) = -h M (T M ) ab <j) M (u) (BIO) 

was used, with M = P, L,T denoting pseudoscalar, longitudinally, and transversely polarized 
vector mesons, respectively. The product h M T M is 



_n-p M 

n M L M — ~ 



if P f lh , (M = P), 

if v % (M — L), (Bll) 
fvVTlaf (M = T). 



B S M 



The coefficients c^ M describe the flavor content of the meson M and are \/2c l 
(—1, —1, 1) for B° s — > (tt°,p°,u), while c^ M = 1 for other decays. 

In order to obtain the corrections from nonperturbative charming penguin to the inclusive 
decay rates, the optical theorem is used in a similar way as in section IIVI To first order in 
a s (2m c ) only the time-ordered product shown in Fig. [Tl 

T c - C =— [ d 4 ze- ip '- z (B\Tf H (z)Q c - c (0)\B) 

m B J 

=4 / d 4 z d% e im ^ l - x ^ z /^B v \T\b v Y n fP L q n ](z) 

J (B12) 

x [(q n )(y)rr M rPLT a (Yic. v -J(0)] [c^ c Y wlfl T a Y^c v - cc ](y) 

x [c v ^ lu P L Y r %](0)\B v ), 

and its hermitian conjugate are needed. The time-ordered product {B\TQ^ tie {z)Q c ^{Qi)\B) 
contributes at order a 2 s {2m c ) and is neglected in our discussion. In Eq. (|B12[) . the time- 
ordered product of the n-collinear fields can be factored out into the jet function 

(0\Tq n (z)Uy)\0) = ip(n.(z-y))5(z ± -y ± ) J ^(-^ j^ {n . k + ie)> (B l 3 ) 
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FIG. 7: Nonperturbative charming penguin contribution to the forward scattering amplitude. The 
blob is the nonperturbative charm contribution and the mirror image is omitted. 



and T c£ becomes 



, dn ■ kdn ■ z 
4 / ■ i 



4tt 



d 4 y ^»{l-XM)-n-k]n-z/2 e in-kn-v/2j m ^ n . J, + ^ 



x (B V \T h Y n (^)^P L ^T M YPLT a Yic. v - c M 



(B14) 



X c_ 



v ^ lfl T a Y^ c (y) ■ c v ,YhvPLYX(0)\B v ). 



If the meson M is a transversely polarized vector meson, Tm = and 7^ vanishes 

because of the spin symmetry. Eq. (jB14j) thus implies that charming penguin effects could 
give a contribution to B — > VrX decays only at sub leading orders in A/m c and/or a s (2m c ). 
We expect that a similar conclusion holds also for the two-body nonleptonic exclusive decays. 
There large transverse polarization fractions, ~ 0.5, have been measured in AS = IB -> 
VV decays (such as B — > 4>K*) that can be charming penguin dominated 42, 3, 0|. This 
may signal substantial l/m c corrections. In naive factorization the transverse component 
on the contrary is expected to be suppressed by 0(rriy /m 2 B ) due to a spin flip. In order to 
explain this large transverse rate,several possibilities of enhanced higher-order contributions 
in l/mfe were suggested jH^, 51 , IH3] - The long-distance charming penguin at leading order 
has also been proposed to contribute to large Rt jl^ . 

For pseudoscalar or longitudinally polarized vector meson, on the other hand, the non- 
perturbative charming penguin contribution is 



r r7 



dn - kdn ■ z 



4tt 



i / d 4 y e i[rn b (l-XM)-n-k]n-z/2 e +in-kn-v/2j m ^ n . fc + ^ 



x (B v \Tb v Y n (?fh^YPLT a Yic- v jO) 



/<u (B15) 

x c^Y^T a Y^Jy) ■ c Vt Y nlu P L YX(0)\B v ). 
The factorization in T cT is more apparent if we rewrite it in a more compact form as 



T- 

* CC 



A 



m b (l-x M ) 



dn ■ I (m,b(l — xm) + n ■ l,n ■ l) J mb (n ■ k + ie) (B16) 
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where we have introduced a new, in general complex, nonperturbative function /, 



(i) 



m h (l-x M ) 



dn-l e in4n - z/2 f£\n-k,n-l) 



2i j d% e m - fc ^ 2 (^|T[(^F n )(^)7^7 1 'P i T a (F 7 Jc_ t „ c )(0 
[c. Vs Y^T a Y^ cc ] (y) [c Vec Y^ lv P L YX] (0)|S„>. 



(B17) 



The integration over n ■ I can be interpreted as the integration over soft fluctuations of b v . 
Taking the discontinuity of the jet function in T c ^ we finally obtain 

-A 



dn ■ I f£ (m b (l — x M ) + n ■ I, n ■ l) 



m b (l-x M ) 



1 



Im J mb (n ■ k + ie) 



7T 



(B18) 



The jet function can be systematically computed in powers of a s (\/Amb). Instead of pursu- 
ing this option, we can treat the convolution of jet function and /jg as a nonperturbative 
function to be determined from experiment. The nonperturbative charming penguin con- 
tribution to the decay rate corresponding to a sum of Fig. [7| and its mirror image is then 



dV^iB -> MX) t . . ,G 2 F f M E 3 M 1 



dE 



T67T a s (2m c 
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N 16vr 2 8r 2 m 2 
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(B19) 



If we include all the possible contributions from the four-quark operators, the nonperturba- 
tive charming penguin contribution to the decay rate at leading order in l/m c ,b and a s (2m c ) 
is written as 

dr- cc (B -> xm) = ql 

8tt 



dE 



M 



-f 2 M m 2 b x 3 M a s {2m c )X^cf M ^ M (l - Ar 2 /x M ) • 2ReT^ ) *^, (B20) 



where the hard coefficients T*f are listed in Table HJ while 
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The nonperturbative function arises from the weak operators with the same color 

, (B22) 



structure as O^-, so that 
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The terms proportional to Air in (jB21j) are smaller than < 2% (< 0.1%) of the terms in the 
first row of (jB21j) for AS = (AS = 1) decays and can be safely neglected. 

The function T C c is independent of the outgoing meson M. In obtaining Eq. ()B20j) an 
expansion in a s (2m c ) was used. If the expansion does not converge one can still parametrize 
the nonperturbative charming penguins by treating the product of a s (2m c ), the LCDA, and 
T C c as a new nonperturbative parameter, to be extracted from experiment. Unlike Tec, 
however, this new parameter depends on M. 
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